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ABSTRACT 


This research investigated multiple detonation diffraction events in order to better 
understand the limits and benefits of diffraction strategies with respect to pulse 
detonation engine design. Hydrogen/air detonations were generated using swept 
ramp obstacles in a 1.27 m long channel with a cross section of 25.4 mm by 88.9 
mm and were diffracted into various multiple-stepped openings. This allowed the 
detonation wave diffraction transmission limits to be determined for hydrogen/air 
mixtures and to better understand reinitiating mechanisms throughout the 
diffraction process. Tests were conducted for area ratios ranging from 2.00-2.60 
with varying equivalence ratios from 0.5-1.5. 

Computational methods were used to better understand the diffraction 
phenomenon using a series of sensitivity studies for different chemistry sets, 
computational cell size and equivalence ratio. Experimental tests used combined 
optical shadowgraph and particle image velocimetry imaging systems to provide 
shock wave detail and velocity information. The images were observed through 
a newly designed explosive proof optical section and split flow detonation 
channel. It was found that area ratios of 2.0 could survive single and double 
diffraction events over a range an equivalence ratio range of 0.8 to 1.14 Area 
ratios of 2.3 survived the primary diffraction event for equivalence ratios near 
stoichiometric for the given step length. Detonation diffraction for area ratios of 
2.6 did not survive the primary diffraction event for any equivalence ratio and 
were unable to transmit to a larger combustor. 
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I. INTRODUCTION 


It is estimated that nearly 80% of the world's energy is derived from the 
combustion of hydrocarbon fuels in various forms. In today's energy conscience 
and profit-driven society, even modest improvements in current systems can 
yield staggering economic gains. The quest to continuously improve energy 
efficiencies of such combustion driven systems has spawned research in 
numerous areas. Pressure gain combustion schemes such as detonation and 
constant volume combustion hold considerable promise in achieving said 
efficiencies [1]. They have inherently higher theoretical thermal efficiencies 
primarily due to their near constant volume detonation based combustion cycle 
which generate much less entropy when compared to constant pressure events, 
creating a pressure rise during heat release. Pulse detonation combustors 
(PDC) exploit such efficiencies. When applied in a pulse detonation engine they 
are estimated to achieve greater specific impulses over many conventional 
systems for a greater range of flight regimes, as seen in Figure 1. 



Figure 1. Comparison of high-speed propulsion technologies (From [2]) 
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Unfortunately, the ability to detonate high mass flow rates in a reliable and 
repeatable manner has been shown to be a tremendous challenge [3]. This is 
especially true for accepted kerosene-based aviation fuels, which are much more 
difficult to detonate than hydrogen. Additionally, the energy required to directly 
initiate a detonation in a combustor large enough to be of practical use in a 
pulsed detonation engine can be very difficult. Naturally, much research has 
attempted to negate complex, high-energy ignition mechanisms by using smaller 
channels or predetonators followed by a transition into a larger section. 
Predetonators or initiators typically involve transitioning a deflagration into a 
detonation through turbulence devices to accelerate the flame. Such devices 
can have losses, so once the detonation is formed, it is normally preferable to 
transition a wave to a larger cross sectional area. This allows for additional fuel 
processing by a surviving detonation wave. Numerous geometries and 
combinations have been used to date, but it remains to be seen whether a 
specific geometry can offer any advantages over others. Hence, the requirement 
to better understand the detonation diffraction phenomena all types of 
geometries. 

Many configurations investigated in the open literature include Brophy et 
al., [4, 5], who studied a number of initiator-combustor combinations, shown in 
Figure 2a through Figure 2c, with oxygen-based initiators with some success. 



Figure 2. Various initiator main combustor conjurations (From [3]) 
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Oxygen-based diffracting initiators were eventually abandoned due to U.S. Navy 
safety concerns and ultimately led the way to a new generation of swept ramp 
initiators, which have been shown to reduce pressure losses by a factor of 5 over 
conventional turbulence devices [6]. 

A great deal of historical detonation research deals with fundamental 
diffraction cases into unconfined spaces and provides tremendous understanding 
of the detonation phenomenon. Unfortunately, there is considerably less 
literature on subcritical reinitiating events, whereby initiator sections are made as 
small possible to support a multiheaded detonation wave front. Such cases are 
of prime concern in the aero propulsion world, since weight and space are at a 
premium. In order to diffract detonations from heavier and more stable fuels with 
larger cell structures such as JP-10, initiators will no doubt be required to 
operate at elevated pressures in a limited amount of volume, so the need to 
understand how aggressive a detonation can diffract is very important. 

This research aims to continue investigations in detonation diffraction by 
exploring multiple diffractions in a 1 to 3.5 aspect ratio channel, in order to better 
understand the limits and benefits of diffraction strategies with this respect to 
pulse detonation engine design. Hydrogen/air detonations were generated using 
swept ramp obstacles in a 1.27 m long channel with a cross section of 25.4 mm 
by 88.9 mm and were diffracted into various multiple-stepped openings. Tests 
were conducted for area ratios ranging from 2.00-2.60, with varying equivalence 
ratios from 0.5-1.5. 

Computational methods were also used to better understand the 
diffraction phenomenon using a series of sensitivity studies for different reduced 
chemistry sets, computational cell size and equivalence ratio. Experimental tests 
used combined optical shadowgraph and particle image velocimetry imaging 
systems to provide shock wave detail and velocity information. 

Results from the research provide baseline design criteria for some 
current pulse detonation engine concepts being investigated in Figure 3. 
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Figure 3. Cross sectional view of a Concentric Tube and multiple high aspect 

ratio channel detonation combustors 


The use of concentric tubes or multiple high aspect ratio channels to transition 
additional fuel to the combustor while isolating the initiator section may allow for 
greater thrust production. In each diffraction plane the detonation front is 
supported by additional air and fuel mixture. Such concepts have the potential to 
produce higher specific impulse values while also smoothing the inherent 
intermittent production of thrust. For the concept to work, the detonation front 
must survive every diffraction event through an appropriate area change. 
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II. BACKGROUND 


A. COMBUSTION WAVE TYPES 

Combustion events can be classified into three categories, explosion, 
deflagration and detonations. Generally, in explosions the rate of heat 
generation is extremely fast and does not require a combustion wave through the 
exploding media. Deflagrations on the other hand propagate at much lower 
velocities with regard to reactants they transmit through. They are essentially 
subsonic expansions waves with pressure drops across the reaction zone due to 
products that are accelerated away from the wave in the opposite direction 
proportional to diffusivity and rate of reaction. Detonations are characterized by 
supersonic combustion waves that propagate extremely fast through a reactive 
mixture increasing the thermodynamic properties sharply, creating a shock- 
induced chemical reaction, which drives the detonation. Some basic detonation 
models include Chapman-Jouguet and the ZND models. Main differences are 
depicted in Table 1. 


Table 1. Detonation and deflagration parameters (From [1]) 



Detonation 

Deflagration 

ul/cl (Mach Number) 

5-10 

0.0001 - 0.03 

u2/u1 

0.4 - 0.7 (deceleration) 

4-6 (acceleration) 

P2/P1 

13-55 (compression) 

« 0.98 (slight expansion) 

T2/T1 

8-21 (heat addition) 

4-16 (heat addition) 

p2/p1 

1.7-2.6 

0.06-0.25 


B. RANKINE HUGONIOT GAS DYANAMIC RELATIONS 

Rankine (1870) and Hugoniot (1887-1889) were the first to investigate 
solutions to the non-reacting conservation equations involving shockwave, often 
termed the Rankine-Hugoniot jump conditions since they applied a control 
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volume across the shock looking solely at the upstream and downstream 
components. In doing so, they were able to derive theses components in terms 
of wave speed, pressures and particle velocity post shockwave. Lee [7] provides 
a nice derivation of these relations, starting with the one-dimensional steady 
state conservation equations across a fixed detonation wave, specified as follows 
with subscripts 0 and 1 representing reactant and product states with velocities 
determined with respect to the wave, as depicted in Figure 4. This theory 
provides a good first approximation to the changes in properties during 
detonations and should also be noted that these equations assume a steady and 
adiabatic front. 






Figure 4. One-dimension steady state model 


Mass itj 

Momentum — pi + 


Energy q + ^ -t- ^ 


Hi 


( 1 ) 

( 2 ) 

( 3 ) 


If we specify an equation of state for both reactant and products using the perfect 
gas law (p = pST), and if we assume with the relationship = if 

and f enthalpy can now be written as follows: 


k = 


r p 

f-ip 


( 4 ) 
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Since we have five unknowns omcfitj), the latter being detonation 

wave speed and only four equations, one more equation is required to close the 
set. 

In order to close this equation, we need to define the Hugoniot curve 
equations. If we start by combining Equations 1 and 2, we obtain: 


Pi -P<i 


_ 2 2 _ 2 2 _ 2 

= 


( 5 ) 



( 6 ) 


In this case, i? = %fp is known as specific volume and rii = pn mass flux per unit 
area. If we start by observing the pressure-specific volume plane of all 
theoretically possible solutions, we can start to better understand relations 
between Rayleigh lines and Hugoniot curves. Looking at Equation 6, we can see 
that real solutions only exist in the detonation and deflagration regions in white in 
Figure 5. 
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/= 


A 

Pn 



Figure 5. Solutions in the p-v plane (From [7]) 


For solutions to exist, we have either >v^ , the case where we 

have a compression solution for a detonation or <p^,-Vii an expansion 

solution and deflagration, as in Figure 5 Lee defines x = = pii/pi and 

7 = doing so one can now rewrite Equation 6 as: 


7 = ( 14 


( 7 ) 



this is the Rayleigh line equation with a slope of we observe that the 

combustion wave velocity is directly proportional to square root of the slope. 
Equation 7 describes a thermodynamic path the gas takes to across the 
combustion wave from state (1, 1) to state (x, y). 
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The Hugoniot curve represents an agglomeration of all possible attainable 
downstream states, for a given specific set of initial conditions, as depicted in 
Figure 6. By eliminating velocities in Equation 3, we obtain the following 
equation for the Hugoniot curve: 


1 

- (^} + IJ) - - {pi ^ P^i) {iTf - 

We also know that solutions to the conservation equations must satisfy both 
Rayleigh and Hugoniot curves simultaneously. In other words, the evolution of 
reactants (x=y=1) to products state (x, y) follow a Rayleigh line until intersection 
on a Hugoniot curve, seen again at Figure 6. For finite heat addition q, the 
Rayleigh line will intersect four points on the Hugoniot curve, two in each of the 
detonation and deflagration branches. In the strong detonation region, both 
pressure and density increase drastically and require overdriven initial conditions. 
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Figure 6. Rayleigh line and Hugoniot curve (From [7]) 


Two other solutions exist when the Rayleigh line is tangent to the 
Hugoniot curve. In this case, we obtain a minimum velocity on the detonation 
branch and a Mach number of the burned gas of Mi=1. The tangency solutions 
are known as Chapman-Jouguet solutions seen in Figure 7. The upper solution 
agrees quite well to the experimental data and is easily calculated using Newton- 
Raphson iterations and has been incorporated in such programs as the NASA 
“Chemical Equilibrium with Applications" (CEA), developed by McBride, Zehe 
and Gordon [9]. 
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C. ZEUDOVICH-VON NEUMANN-DORING (ZND) ONE DIMENSIONAL 

WAVE STRUCTURE 

A simple and widely accepted model of a detonation wave is the ZND 
model comprising of a strong shock front followed by a much thicker induction 
and reaction zone travelling at Chapman-Jouguet velocity. The shockwave 
initially adiabatically compresses the mixture to a peak pressure known as the 
Van Neumann state, enabling dissociation of the molecules within a thermally 
neutral induction zone. At a certain point, enough free radicals are created to 
start the chain branching reactions eventually converting all reactants into 
products within a reaction zone. The rapid chemical reaction rates enable a very 
quick rise in temperature dropping both pressure and density in the reaction zone 
the resulting expansion accelerate the gases away behind the front coupling it 
the acoustic front [7]. 
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Figure 8. Model and Real detonation (From [8]) 

D. DETONATION CELL STRUCTURE 

Detonation fronts in real life have complex three dimensional structures. 
The wave front contains multiple and constantly interacting curved incident, Mach 
stems and reflected shock waves meet at so called triple points. These 
structures are highly transient due to spatial and temporal instabilities. If we 
observe the highly resolved Schlieren photograph on Figure 9, we can see what 
the reaction looks like before two transverse waves collide. The accompanying 
sketch provides some additional explanation into all the interconnected events. 
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Figure 9. Schlieren photograph of a detonation in a thin channel (From [7]) 

The shock front has multiple colliding waves that progress transversely 
across the front, creating a triple point path that forms a diamond shape. The 
chemical reactions are not always perfect during a detonation event; note the 
turbulent nature of the reaction zone behind the detonation with numerous shear 
layers and unburned pockets that are consumed. In highly unstable mixtures like 
hydrogen-air and many hydrocarbons have strong transverse waves which tend 
to drive the detonation. This can reduce the amount of chemical energy 
contributing to the main front leaving unburned pockets to be consumed later by 
turbulent diffusion. These mixtures can also exhibit irregular cell structures and 
irregular triple point trajectories. Contrary to some highly regular cell structured 
mixtures that contain heavy amounts of argon as depicted in Figure 10. 
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Detonations with regular cell structure are less reliant on transverse wave 
collisions to continue the propagation of the front. Generally, these mixtures are 
more reactive. 



Triple point trajectories 


Mach stem 

Chemical induction layer 
Incident shock 


Shear layer 


Transverse shock 


Figure 10. Example of a regular triple point cell structure of a detonation 

(From [7]) 


Triple point trajectories can be mapped by soot foil techniques whereby 
shear and pressure gradients rearrange soot particles leaving behind evidence of 
triple point paths. From these soot foils traces, cell width A can be determined 
and provides for an estimation of mixture reactivity. The smaller the cell width is 
the tighter the coupling is between chemical kinetics and gas dynamic effects. It 
should be noted however that there is a large variance in measurements in cell 
structure via soot foils as observed by Shepherd [10]. Soot foil measurements 
can be very subjective with errors up to 50%. 
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E. 


STOICHIOMETRY 


Combustion conditions of fuel-oxidizer mixtures are generally specified in 
terms of equivalence ratio O. This characterizes the proportion of fuel to oxidizer 
in the mixture. Defined as the ratio of the fuel-oxidizer ratio to the stoichiometric 
fuel-oxidizer ratio and can be mathematically represented as follows: 


4 > 


^ Juel oxidizer 




idizer ) 


^ fiiel j ^oxidizer 


stoichiomenic 


Juel I ^oxidizer ), 


sioichionettic 


( 9 ) 


The composition and thermodynamic properties of fuel-oxidizer mixtures during 
combustion and detonation are dependent on the equivalence ratio. 
Consequently the detonation cell width is also highly dependent on mixture fuel 
ratio. A ratio greater than 1 implies excess fuel or fuel rich mixtures than would 
be required for a stoichiometric reaction and conversely a ratio less than one 
implies excess oxidizer or fuel lean mixtures. Looking more closely at the 
dependence of detonation cell width by inspecting Figure 11, we can see a very 
strong relation to equivalence ratio. Cell width for a particular mixture is smallest 
at or very close to 0=1. 
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Figure 11. Variation in detonation cell size with fuel air composition (From [1]) 

F. DETONATION DIFFRACTION 

Diffraction occurs when a confined detonation is allowed to expand into a 
larger cross sectional area at which point expansion waves are communicated to 
the flow from the behind the front. Diverging transverse waves affected by this 
perturbation produce unsteadiness and curvature in the front, which by definition 
also reduces chemical reaction rates. This reduces the energy is release rate 
thereby providing less support to the detonation front. Lee illustrates this in 
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Figure 12, we see a rarefaction fan eventually propagates toward the axis of the 
channel or tube and the flame front decouples from the acoustic front in a total 
failure. 



Figure 12. Rarefaction fan penetrating propagating a toward detonation axis 

(From [8]) 


Whether a detonation survives or fails an expansion depends on a number of 
factors: mixture composition, thermodynamic initial conditions, confining 
geometry or boundary conditions and detonation velocity at diffraction plane. 
Three different cases have been observed: super critical where a detonation 
successfully transitions a diffraction plane, subcritical cases occur when a 
detonation wave is quenched and critical cases are known as the transition point 
between the two previous cases, seen in Figure 13. 
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(a) Supercritical: suc¬ 
cessful traiismissioii. 


(b) Critical; re-initiation 
as shock decouples 
from reaction zone. 


(c) Snbcritical: complete 
failure as shock decou¬ 
ples fioiii reaction zone. 


Figure 13. Diffraction regimes (From [8]) 


The detonation wave diffraction phenomena is not new; however, the last 
75-plus years of research in this area has still has not yielded a complete 
understanding of all the fluid dynamic complexities in this event. With all 
detonations, high pressures and temperatures, short and varying time scales 
have stymied much of the research involved in this area. Lafitte [11] was the first 
to attempt the transition of a detonation wave from a 7mm tube into a spherical 
chamber using a mixture of CS 2 +O 2 . He was unsuccessful. Zeldovich's [12] 
Streak camera experiments exhibited detonations that decayed completely into 
flames for some conditions while others continued to detonate across the 
diffraction plane, he established that for various tube diameters a critical diameter 
exists for successful reinitiation into a larger space, based on a mixtures initial 
conditions and confinement properties. With improvements in experimental 


18 




techniques in the 1950s through 1970s, cellular structure of the detonation wave 
was better characterized. Shchelkin and Troshin [13], were the first to relate cell 
size to the Zeldovich-Newmann-Doring (ZND) model reaction length. 
Mitrofanov and Solouhkin [14] noted the first empirical correlation relating 
detonation cell width to critical tube diameter; they proposed that the critical 
diameter is 13 times the detonation cell width A. They demonstrated that for a 
number of hydrocarbon fuels that if a detonation front had 13 or more cell widths 
across a tube diameter, the detonation would survive a diffraction event into any 
open space with the same chemistry. Follow-on research [15, 16, 17] finds this 
theory breaks down for more regular cell structured fuel mixtures and higher 
aspect ratio rectangle configurations as seen in Figure 14. Extreme aspect ratio 
channels great than 5 asymptote to critical channel widths Wc of 4 times the 
detonation cell width A. Such channels appear to more easily transition 
detonation waves as other transverse modes are suppressed and expansion is 
restricted to two dimensions. 



Figure 14. Critical channel height as a function of aspect ratio (From [7]) 
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Cell structure and chemical kinetics, as well as initial and boundary condition, 
play a major role in how detonations waves expand, break down and reinitiate, 
as well as the initial and boundary conditions [17, 18]. 
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III. COMPUTATIONAL SETUP 


A. INTRODUCTION 

Computational fluid dynamic (CFD) software was used extensively during 
the research in order to compliment experimental results and to better 
understand all the complexities of a diffracting gaseous detonation wave. 
Numerical results also provided additional details of the flow field that were 
difficult or impossible to do so experimentally. Simulations were largely based on 
solving the two dimensional non-linear compressible inviscid reactive 
conservation equations known as the Reactive Euler equations below in 
conjunction with multi-step reduced chemistry sets with nine species (H 2 , O 2 , H, 
O, OH, HO 2 , H 2 O 2 , H 2 O, and N 2 ) and 18/19 elementary reactions in rectangular 
coordinates.. 
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In short conservation of mass, momentum in both x and y directions, energy 
along with all species equations in terms of Arrhenius rates and an equation of 
state were used. Metacomp’s CFD++ software was used to solve this equation 
set. Structured meshes were generated using Pointwise meshing software 
package. In order to set-up a baseline for the study, a number of sensitivity 
studies were conducted using computational cell size, degree of overdrive, 
stoichiometry and chemistry. Sensitivity simulations were obtained using both 
implicit and explicit schemes. Due to the transient nature of detonations, implicit 
simulations, a CFL of 1 was used to ensure maximum stability. In all explicit 
time-marching schemes, CFL number was ramped from 0.5-1 during the first 
1500 iterations allowing the simulation to resolve initial deflagration to detonation 
formation. Consequently, the time step selection in explicit time-marching 
simulations were critical for stability, in our case as a detonation wave crosses a 
discrete grid, the time step must be small enough to resolve it. 

B. REDUCED CHEMISTRY 

In order to correctly predict where and how detonation waves would 
diffract and whether or not the end state survived or failed, heavy reliance on 
reduced chemistry was required. Simple H 2 and O 2 reactions systems have 
been known to encompass hundreds of reactions. However, computational 
resources are normally limited. There has been significant research in using 
reduced chemical models that capture the most important characteristics with 
respect to a particular application. Models in the past have either assumed one- 
step equilibrium based chemistry having limited success in predicting 
detonations. This has tended to over simplify the reaction and over predict 
product temperature. Computational power today has allowed for multistep finite 
rate chemistry allowing better prediction of heat release and ignition delay times 
(induction time) during a detonation. 

This research looked at three reduced sets that have less than 20 
reactions, none of which were developed exclusively for detonation based 
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combustion. The first set developed by Kurkov & Burrows at NASA for scramjets 
was based on lower combustion pressure regimes [19]. The second set was a 
modified Kurkov set that included the some energy absorbing NO chemistry 
based on a reduction scheme from Kuo [1]. The last set used was developed at 
Stanford by Peterson and Hanson for ram accelerators for a range of higher 
pressures for up to pressured of 50 atm [20], as seen in Table 2. Hydrogen air 
detonations as explained earlier are highly unstable reactions and have 
tremendous pressure fluctuations that can range from 10 atm to 27 atm or from 
0.5Pcj to 1.5Pcj. [7]. The reaction rates were expressed in Arrhenius form as 
follows: 



(18) 


The reaction rate constant, k was a function of temperature, pre-exponential 
parameter A described the frequency of collisions for the species that result in a 
reaction as well as the preferred orientation of the molecular species that favor a 
reaction. The exponent n expresses the temperature dependence of the pre¬ 
exponential parameter. Activation energy, Ea defines the minimum amount of 
energy required for a reaction to occur and R is the universal gas constant. 
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Table 2. Peterson & Hanson reduced set for high pressure H 2 /O 2 , N 2 
reaction mechanism for ram accelerator (From [20]) 


Rate coefficient*’ 


Number 

Reaction® 

A 

n 

E 

CojiiJiients^ 

1 

0 + H2 = H + 0H 

5.00 X Kf 

2.70 

6,290 

(R3) 

2 

H + 02 + M = H02 + M 

2.80 X 10 '® 

-0.90 

0 

(RSS)** 

3 

H + O 2 O 2 = HO 2 + O 2 

3.00 X 10^® 

-1.70 

0 

(R34) 

4 

H + O 2 + H 2 O = HO 2 + H 2 O 

9.38 X 10'® 

-0.80 

0 

(R35) 

5 

H + O 2 + N 2 = HO 2 + N 2 

2.60 X 10'® 

-1.20 

0 

(R36) 

6 

H + 02 = 0 + 0H 

8.30 X 10'3 

0.00 

14,413 

(R38) 

7 

H + HO 2 = O 2 + Hp 

2.80 X 10'3 

0.00 

1,068 

(R45) 

8 

H + HO 2 = OH -h OH 

1.34 X lO'-' 

0.00 

635 

(R46) 

9 

H + H 2 O 2 = HO 2 -h H 2 

1.21 X 10® 

2.00 

5,200 

(R47) 

10 

OH -h Hp = H 2 O + H 

2.16 X 10® 

1.50 

3,430 

(R84) 

11 

OH -h OH + M = H 2 O 2 + M 

7.40 X 10'® 

-0.40 

0 

k,„f, (R85^f 



2.30 X 10'® 

-0.90 

-1,700 

Ao 

12 

0H-hH02 = 02 + H20 

2.90 X 10'® 

0.00 

-500 

(R87) 

13 

OH -h H 2 O 2 = HO 2 + H 2 O 

1.75 X 10 '® 

0.00 

320 

A„,tR88)E 



5.80 X lO'"' 

0.00 

9,560 

Ai„(R89)s 

14 

HO 2 + HO 2 = O 2 + H 2 O 2 

1.30x 10" 

0.00 

-1,630 

Ar^.tRllS)'’ 



4.20 X 10'^ 

0.00 

12,000 

Arf, (R116)'’ 

15 

0 + 0 + A/=02+it7 

1.20X 10'® 

-1.00 

0 

(Rl)’ 

16 

0 + H + A/=OH + M 

5.00 X 10'® 

-1.00 

0 

(R2)" 

17 

H + 0H + A/ = H20+A/ 

2.20 X 10®® 

-2.00 

0 

(R43)i 

18 

H + H + M = H 2 + M 

1.00 X 10'® 

-1.00 

0 

(R39)'‘ 


Note; species = H* O, OH, H 2 , O 2 , HOi, H 2 O 2 , H 2 O, and N 2 (or Ar, He, etc.X 
“All reactions are reversible. 

= AJ"exp(—£/ RTy, units are in cal, mol, cm^, and s. 

Corresponding RAMEC^*^ reaction numbers are in parentheses. 

M does not include O 2 , H 2 O, or N 2 ; all collision efficiencies = 1.0. 

'Collision efficiencies for iW; N 2 = 1.0, H 2 = 2.0, H 2 O = 6.0, and Ar = 0.70,all others = 1.0. 
*^See Table 2 for Troe parameters. 

^Rate coefficient is non-Arrhenius; A -13 = A„ -l-Afr. 

•'Rate coefficient is non-Arrhenius;A |4 = -l-Aj. 

‘Collision efficiencies for Af: N 2 = 1.0, H 2 = 2.4, H 2 O = 15.4, and Ar = 0.83, all others = 1.0. 
^Collision efficiencies for M; N 2 = 1 . 0 ,H 2 = 0.73, H 2 O = 3.65, and Ar = 0.38, all others = 1.0. 
Co Hi Sion efficiencies for jW; N 2 = 1.0, H 2 = 1.7, H 2 O = 7.0, and Ar = 0.63, all others = 1.0. 


C. COMPUTATIONAL MESH 

The computational mesh used was a two dimensional channel with an 
open step in the positive y direction with a 1:6 height expansion ratio. A 
resolution of 50 pm or roughly 300 computation domains in one detonation cell 
for H 2 was used equaling to approximately 4.5 million cells. The mesh was 
developed using Pointwise 64-bit meshing software and converted to run in a 
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CFD++ environment, as seen in Figure 15. Some researchers have reported 
that the resolution requirement to fully resolve complex chemistry simulations is 
an element size of only 0.1 pm for the CJ H2-air detonations at 1 atm and 298 K 
[21]. Essentially three orders of magnitude finer of the induction zone thickness. 
Such resolutions would be well beyond currently available computational 
resources and could not be attempted in this study. Additionally, previous 
studies at NFS [22] have shown that this resolution predicts detonation velocity 
within 5-10% of CJ values. It was hoped that this resolution would be sufficient 
to capture failure modes of a detonation diffracting around an open step. 

The model was scaled down slightly from the original experimental setup, 
with channel height scaled down to 1:2 and the channel length scaled down 1:5, 
this allowed for solutions within a week. 
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Figure 15. Computational domain (in m) for diffraction studies 













D. 


BOUNDARY AND INITIAL CONDITIONS 



Figure 16. Boundary and Initial Conditions Computational results 


Detonation waves in the simulations were created by introducing 
overdriven high pressure products in small region at the beginning into a main 
body of reactants. Detonation products were initiated at two times CJ Pressure 
and at CJ Temperature conditions, which as alluded to above was found to be 
the minimum overdrive found for given conditions to initiate a detonation wave 
characteristic of the fuel-oxidizer mixture as part of the initial conditions. Mass 
fraction of the detonation products were obtained from NASA Chemical 
Equilibrium Analysis (CEA) [9] and a sample is included in Appendix A. Reactant 
initial conditions were assumed to be at 1 atm and 298 K, essentially matching 
experimental conditions for further validation. Typical run set structure was as 
follows: 
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Equation Set Type 

• Unsteady Compressible Euler 

• Equation of State: Ideal Gas 


Riemann Solver 

• Minimum Dissipation: LHS 

• No pressure switching 

Boundary Conditions (See Figure 16) 

Time Integration 

• Implicit CFL=1 

• explicit (for some runs) based a CFL ramped from 0.5 - 1 
over the first 1500 iterations 

Spatial Discretization 

• 2nd Order Accuracy in Space 

• Dimension of polynomial: 2-D XY 

• Axis stability enhancement: Yes 

• Types of Total Variation Diminishing (TVD) Limiter: Minmod 
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IV. EXPERIMENTAL SETUP 


A. DETONATION CHANNEL MAIN COMPONENTS 

The combustor design was driven by three main requirements, the desire 
to use hydrogen and air as reactants, a large field of view in the test section in 
order to facilitate both shadowgraph and particle image velocimetry and lastly, 
the ability to adapt the existing main tube design in test cell one. The main 
components of the detonation diffraction test setup consists of an inlet manifold, 
a 48-inch long rectangular channel in which an inner split flow channel is inserted 
with a one 3.5-inch tall by 1-inch wide inner channel dimensions and an optical 
section. This configuration allows for detonations to pass through a lower section 
and diffract through an upper section carrying unreacted fuel-air mixture into a 
viewable optical section seen at Figure 17. Appendix B includes all engineering 
drawings for the split flow insert and optical section components as well as the 
new manifold 



29 


























































































Figure 18. Main test section components 


B. COMBUSTOR OPERATION 

This pulse detonation combustor setup operates on a relatively simple 
single cycle mode in order to facilitate optical experimental observations with 
shadowgraph and Particle Image Velocimetry (PIV). It was operated on an 
optics table consisting of two electro pneumatic ball valves that controlled fuel 
and seeding air along with a main air duct regulated by a computer controlled 
pressure regulator. Also, various check valves, two hand pressure regulators to 
set precise pressures for fuel and seeding along with two mass metering chokes. 
Two lasers were employed to illuminate the detonation front for image capture for 
their respective camera system using a continuous Argon-Ion for shadowgraph 
and a double pulsed Nd-Yag for PIV, seen in Figure 19. Two kistler high-speed 
pressure transducers acquired high speed data and were used to determine 
detonation velocity as well as trigger the PIV system via a BNC 555 pulse delay 
generator. Labview commercial control and automation software was used to 
operate the test engine. The standard sequence consisted of filling the upper on 
lower channel with a seeded titanium dioxide hydrogen/air mixture and igniting it 
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and subsequently measuring pressure, detonation velocity, acquiring high speed 
shadowgraph images and particle image velocimetry images. See Appendix C 
for standard operating procedures. 



Figure 19. Overall experimental setup 
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Figure 20. Overall experimental Setup 


C. OPTICAL SECTION 

The tube was attached to an 8.5-inch tall by 12.5-inch long and 4.5-inch 
wide optical section with opposing 5-inch by 9-inch viewing areas using 45- 
degree chamfered flush fit BK7 optical glass. Various step configurations can be 
installed into the section. The design was based on a previous section by 
Fludovich [4], however the new section was designed to be more modular with a 
larger optical viewing area bolting directly to the detonation channel with access 
ports above and below, seen in Figure 21. The entire section was machined into 
mainly one inch thick 304 stainless steel sections and was assembled with 
socket cap head screws. The strong design is required in order to withstand any 
bending that could compromise the 1.75 inch BK7 optical glass it was designed 
to encase. 
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Figure 21. Stainless Steel 304 optical section 


Finite element analysis was performed on the BK7 glass components of 
the older and newer designs to ascertain maximum stress and deflection under 
shock loading conditions in order to ensure that the new design would perform 
equally or better. A glass thickness of 4.45 cm (1.75 in) was selected to 
accommodate the detonation impulse pressures of up to 36 atm for up to 1 ps 
shock loading. Results of these analyses are presented in Figures 22 and 23. 
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Figure 22. Deflection results for BK7 glass 
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Figure 23. Principal stress results for the BK7 glass 
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D. 


LOWER CHANNEL SPLIT FLOW INSERT 


In order for a detonation to diffract into an open step using the existing 
main square channel, a split flow channel assembly was designed to easily slide 
into the existing hardware, essentially separating the flow into two channels, an 
upper and lower allow for a detonation to initiate. The detonations were initiated 
in the lower section expanding into an adjacent top channel at the diffraction 
plane directly observed through the optical section. The flame was transitioned 
into a detonation using 10 interspaced 30-degree swept ramps with a total 
blockage ratio of 40%. This section was machined from 6061 Al stock and is 
socket cap screw assembled, more detailed drawings can be viewed in 
Appendix B. 



Figure 24. Split flow insert. 
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E. 


IGNITION SYSTEM 


The detonation is initiated by a Unison Vision-50 variable capacitative 
discharge using an aviation grade spark plug mounted at the head end of the 
tube. Prior to testing, the igniter was configured manually to provide a 2 Joule 
spark and was remotely triggered by the Labview software (Figure 25). 



Figure 25. Unison variable capacitative system (From [23]) 

F. FUEL AND AIR DELIVERY SYSTEM 

High pressure air was provided to the test cell from the rocket lab facility 

air system with hydrogen provided from a single bottle. Main control room 

computers controlled supply pressures for both air and hydrogen with the 

Tescom ER3000 Version 2.0 software, allowing reactant pressures to be set 

remotely as required. Fuel, air and seed delivery revolved around the main six 

port inlet manifold which was designed to adapt to the existing pair of 1 inch main 

air lines. Further upstream air is fed through a 1.5-inch pipe using a 0.235-inch 

flow metered choke. The manifold design delivers a high enough mass flow, 
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ensuring swept ramp effectiveness and thus reliable detonations. The design 
allows for low pressures seed particulate penetration into the main flow and 
ensured the seed reservoir's 2 bar pressure limit was not exceeded. Main air 
lines were broken down into two successive pairs of half inch lines where fuel 
and seed was fed to the main flow in two stages into the six-holed manifold that 
accepted up to five lines with one hole tapped for a standard turbine plug., as 
depicted in Figure 26. Hydrogen was supplied using half-inch tubing that in turn 
fed two main air arms just prior to the four line split. 



Figure 26. Inlet manifold and feed lines 


Unfortunately hydrogen and seeding could not be regulated using the 
ER3000 systems since operating pressures were lower then what this system 
was capable of delivering to within reasonable errors. Hand ball valve regulators 
were used instead, since typical operating pressures were under 100 psig and 10 
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psig, respectively. All pneumatic ball valves were controlled by LabVIEW 
software via Crydom control solenoid switches located in an electronics cabinet 
located in the test cell as seen in Figures 27 and 28. Choke diameters used 
during this testing were 0.236 inch for the air supply and 0.157 inch for the 
hydrogen supply. Two check valves located in the fuel and seed lines were 
installed to prevent the possibility of backflow due to the detonation event. 
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Figure 27. Fuel and air delivery system (From [23]) 


G. INSTRUMENTATION 

Instrumentation employed all through the testing controlled initiation as 
well as data collection. Each of the pressure signals were routed through one of 
four National Instruments (Nl) 14-bit PXI-6115 cards mounted in the Nl PXI- 
1000B chassis, shown in the upper half of Figure 28. This chassis interfaced 
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with the computers in the control room through the Nl PXI-MXI-4 PXI Bridge and 
was capable of collecting either real-time or high-speed buffered data. Optical 
data imaged from the high-speed camera and PIV systems were routed to a 
dedicated desktop computer in the control room. 



Figure 28. National Instruments PXI-1000B Chassis (upper section) and 
Crydom Control Solenoid Switches (lower section) (From [23]). 

H. DYNAMIC PRESSURE TRANSDUCERS 

Inside the channel, two high-frequency Kistler dynamic pressure 

transducers were mounted as a pair in order to calculate the wave velocity based 

on the pressure transients. The transducers were connected to Kistler Type 

5010 Dual Mode amplifiers which route the data to two Nl PXI-6115 data cards 

in the PXI-1000B for high-speed data collection. Once testing was initiated, the 

cards were configured to begin collecting data with the triggering of the 
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capacitive discharge igniter and collected data at a rate of 500 kHz. Data was 
stored via card buffer and later saved to computers located in the control room. 
Pressure data supplied from Kistler transducers provided high-resolution 
confirmation of shock passage. 

I. CONTROL 

Every test sequence was initiated and controlled using Labview software. 
All pressure transducers were linked into to the code to allow for active 
monitoring during every test sequence. A sample control panel can be seen in 
Figure 29. 


Test Cell 1 Live 

Filename and path 
I C:\JuilletData\ldecl0.txt 
Highspeed file path 
I C:\JuilletData\HighSpeedldecl0.txt 





LmJ 




Kistlerl 

Kistler2 




Sample Rate (hz) 
1000 

Samples 

10 

Min value (volts) 
0 

Max value (volts) 


H2 Fuel 
Pressure (psig) 


Main Air PreChoke 
Pressure (psig) 


Node4 Seeder Air 
Pressure (psig) 


Sample High Speed Rate (hz) 
500000 

Samples Highspeed 
200000 

Min value HS (volts) 

-10 

Max value HS (volts) 

10 


Figure 29. Sample Labview control panel 


40 








J. SHADOWGRAPH 

Opposing optical windows allow for the transmission of a shadowgraph 
light source to pass through the section during a detonation. Shadowgraph setup 
utilizes an Ion Laser Technology Model 55001 Argon Ion laser emitting 750 mW 
beam at 488 nm passing into a 6.3 pm spatial filter placed in front of the laser 
aperture. This produced a spatially filtered light source, which became 
collimated after being reflected from a concave mirror. The resulting beam was 
directed through the optical test section using a flat mirror. On the opposite side 
of the optical test section, a flat mirror directs laser light toward a focusing mirror 
placed which in turn is directs the light directly into the ultra 17 intensified high¬ 
speed camera. 

The main rationale for choosing shadowgraph as opposed to Schleiren 
technique was due to the high coherence length of the argon ion laser system. 
Such high coherence length, do not provide for good Schleiren images. The 
argon ion system was used as monochromatic light source. 

Invisible Vision's Ultra 17 high-speed camera suite was used to capture all 
shadowgraph images. The system was triggered from the PI kistler high speed 
pressure probe with delays ranging from of 225-275 ps depending on what 
portion of the optical window was being captured. At the 150,000 fps mode, a 
resolution of 512 X 512 pixels was possible with this system. Exposure settings 
varied between 50-500 ns, depending on the amount of light in the test cell while 
allowing gain to be lower than 30%. Inter frame times were 6.17 ps. 
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Figure 30. Shadowgraph setup 


K. PARTICLE IMAGE VELOCIMETRY (PIV) 

The PIV arrangement used a Litron Nd:YAG Dual Cavity pulsed laser 2 x 
200 mJ at 532 nm with up to 15 Hz pulse rate with variable light optics. The laser 
sheet was passed underneath perpendicular to the flow through a narrow optical 
slot Image as seen in Appendix B. Image capturing was done using a Lavsion 
Imager Pro X 2M 1648 x 1214 pixel resolution camera fitted with a 50 mm f1.4 
Nikon lens and 532 nm filter. Laser and camera timing was accomplished using 
a PC based programmable timing unit (PTU-9). All image post-processing was 
completed with DaVis software installed on a standard quad core PC. 

In order to properly seed, a reacting and detonating proper seed material 
would be required to withstand high temperature ad pressures while also 
exhibiting minimum interference to the flow. Micron particle Titanium Dioxide 
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was selected for its good stability at high temperatures and its relative inertness. 
As eluded to earlier manifold flow seeding was accomplished using the Lavision’s 
Particle blaster 100 and titanium dioxide seeding particles. Since testing 
involved simultaneous PIV and shadowgraph, the PIV camera was angled out of 
the shadowgraph optical path, requiring calibration within the DaVis software 



Laser sheet 


Double Frame 
PIV camera 




Double Pulse 
ND;Yag laser 


O 


Figure 31. PIV setup 


L. OPTICAL SECTION TEST CONFIGURATIONS 

There were five major test configurations: Configuration 1 consisted of 
open diffraction with no step sections and a series of cold flow tests. 
Configuration 2, 3 and 4 were single step diffraction setups of various height 
ratios; 5 and 6 were double step diffraction setups with two different height ratios. 
Figure 32 depicts common setup configurations. 
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Figure 32. Single and double step configurations 
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V. COMPUTATIONAL RESULTS 


A. OVERVIEW OF RESULTS 

It was originally envisioned that the CFD portion of this thesis could 
potentially shed light and perhaps help us better understand the complexities of 
gaseous diffraction detonation wave phenomena while also assisting in the 
design of some of the components. Simulations in general were effective in 
forming a detonation given the various parameters but once the detonation front 
was formed, failure of the front during diffraction was not attained. Unfortunately, 
not one simulation was observed to completely fail across the diffraction plane 
and thus precluded us from any useful results in computational diffraction 
detonation. 

Two main reasons were suspected for simulations that could not resolve a 
failing detonation front, namely poor mesh resolution and the reduced chemistry 
sets utilized. Poor Resolution (causing artificial reactivity) resulted in an 
unresolved the triple shock system at the front. Incorrect chemistry could not 
predict the proper heat release and ignition delay time (induction time), as well as 
not enough diffusive elements in the simulation. This resulted in solutions that 
were overly reactive and not capable of failure during a diffraction event. 

Simulations for the smallest mesh size of 0.05 mm used for the study ran 
for roughly 240 hours on 128 processors on the new NPS hamming Linux 
cluster, while the largest mesh sizes took only 24 hours to run with the same 
computational settings. CFL number for all simulations was in between 0.5 and 
1 and with varying time steps depending on the mesh size. Table 3 is a 
summary of all the computed cases. 
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Table 3. 


Summary of computed results with ATM initial conditions 


Run # 

Mes 

h 

Size 

(mm) 

4> 

Reduced 
Chemistry Set 

Numericai Scheme & other 
settings 

Result 

1 -23jul 

0.05 

1 

Kurkov Burrows 

Implicit, CFL 1 

No failure 

2 -26jul 

0.1 

1 

Kurkov Burrows 

Implicit, CFL 1 

No failure 

3- 

16aug 

0.1 

1 

Kurkov Burrows 

Explicit, CFL ramped 0.5-1 

No failure 

4- 

20aug 

0.05 

1 

Kurkov Burrows 

Explicit, CFL 1, 1C 1.5Pcj 

No 

detonation 

5- 

23aug 

0.05 

1 

Kurkov Burrows 

Explicit, CFL 1, 1C 2Pcj 

No failure 

5- 

23aug 

0.05 

1 

Kurkov Burrows 

Explicit, CFL 1, 1C 2Pcj 

No failure 

6- 

27aug 

0.05 

1 

Kurkov Burrows 

Implicit, CFL 1, LHS 

dissipation 

No failure 

7- 

30aug 

0.05 

1 

Kurkov Burrows 

Implicit, CFL 1, LHS 
dissipation, 1C 1.5Pcj, 

viscosity on walls 

No 

detonation 

8- 

31aug 

0.05 

1 

Kurkov, Burrows 

Implicit, CFL 1, LHS 
dissipation, 1C 2Pcj, viscosity 
on walls 

No failure 

9-2sep 

0.05 

1 

Modified Kurkov 
Burrows 

Implicit, CFL 1, LHS 

dissipation, 1C 2Pcj, added 
NO chemistry 

No failure 

10- 

7sep 

0.05 

1 

Modified Kurkov 
Burrows 

Implicit, CFL 1, LHS 
dissipation, 1C 2.5Pcj, 

viscosity on walls, added NO 
chemistry 

No failure 

11- 

lOsep 

0.05 

0.8 

Modified Kurkov 
Burrows 

Implicit, CFL 1, LHS 
dissipation, 1C 2Pcj, added 
NO chemistry 

No failure 

12- 

loct 

0.05 

0.8 

Modified Kurkov 
Burrows 

Implicit, CFL 1, LHS 

dissipation, 1C 2Pcj, added 
NO chemistry, symmetry 
plane 

No failure 

13- 

2oct 

0.05 

0.6 

Modified Kurkov 
Burrows 

Implicit, CFL 1, LHS 

dissipation, 1C 2Pcj, added 
NO chemistry, symmetry 
plane 

No failure 

14- 

12oct 

0.05 

0.6 

Modified Kurkov 
Burrows 

Implicit, CFL 1, LHS 

dissipation, 1C 2Pci, added 

No failure 
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NO chemistry 

15- 

0.05 

1 

Modified Kurkov 

Explicit, CFL ramped from 0.5 

No failure 

26oct 



Burrows 

-0.95, ran faster, 2Pcj, added 
NO chemistry 


16- 

0.05 

1 

Hanson, 

Explicit, CFL ramped from 0.5 

No failure 

5nov 



Peterson 

-0.95, ran faster, 2Pcj, added 
NO chemistry 


17- 

16nov 

0.05 

1 

Hanson, 

Peterson 

Explicit, CFL ramped from 0.5 
-0.95, ran faster, 2Pcj, added 
NO chemistry, symmetry 
plane 

No failure 


B. DETONATION FRONT STRUCTURE 

The following figures show the computed contours of pressure, OH mass 
fraction, pressure and temperature plots for 0.1 mm and 0.05 mm mesh sizes. 
The combustion front at these resolutions did not capture enough front details. 
Best resolved fronts came from the OH contours at the finest resolution. 

Upon first inspection at the detonation front on the OH mass fraction 
contours right before the diffraction plane, the front appears to lack real 
sharpness of a triple point that we would normally from a two dimensional 
simulation. Very faint keystone structures were only observed with a 0.05 mm 
mesh. There was no defined zone where the three shocks meet. Looking at 
Pintgen’s extensive PHD work using planar laser induced florescence (PLIF) of 
the OH radicals he revealed their presence near the front in the distinct keystone 
patterns is clear in Figure 33. 
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Figure 33. Observations in the sub-critical regime, combined Shlieren, OH* 
and PLIF imaging of the highly irregular H 2 -N 2 O mixture (From [18]) 


The H 2 N 2 O mixture in Pintgen work has a much more irregular cell 
structure and is much more reactive than our H 2 air case with cell sizes 
averaging 1-4 mm. We use the case only to compare irregular and subcritical 
type trends. Normally during high temperature and pressure events, dissociation 
occurs immediately after shock passage resulting in high mass fractions in 
vicinity of the front, this is just barely apparent with these contours. A run with a 
5pm mesh size or smaller would have been attempted except that but such runs 
would have taken months for results with NPS' current computational resources 
and would not have been feasible given current time constraints. It is understood 
that for structured meshes a 5 pm resolution has been reported as the starting 
point for adequate detonation front structure [21]. Looking at Figure 34 OH mass 
fraction contour plots, there are drastic improvements in front details moving from 
100 pm to a 50 pm structured mesh—regrettably, not enough to resolve any 
jagged edges or keystones that we can clearly see in Pintgen's results. 
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Figure 34. Mass fraction OH contour plots of detonation front using a 100 pm 

and 50 pm structured mesh 


With the pressure contours for all three mesh sizes, we can measure 
detonation cell size A from remnant transverse pressure waves right behind the 
front. It is interesting to note that average cell size in each mesh size seem to be 
consistently smaller then experimental values of 1.5 cm for Hydrogen. For our 
baseline case with a mesh size of 50 pm, we found an average cell size of 
0.22 cm, indicating a much more reactive mixture. Many researchers have 
reported that in order to computationally acquire accurate cell size, proper grid 
resolution is critical. Poor resolution has also been known to cause smaller than 
normal detonation cell sizes and in other words a much more reactive mixture. 
Given the computational constraints we could not confirm this since we did not 
run the 5pm mesh case. Transverse waves also revealed some additional 
irregularities; they seemed to persist much longer after the front had passed than 
normal. They occurred much more regularly then one would expect for such an 
irregular mixture like H 2 -Air. Pressures in and around the front were on average 
20% higher for both Von Neumann and CJ pressures. This was expected since 
nitrogen chemistry was not included. 
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Figure 35. Pressure contour plots of detonation front using a 100 pm and 

50 pm structured mesh 


Synonymous with detonations was also the turbulent nature of the 
reaction zone. Viscosity and transport effects played an extremely vital function 
within the small scale phenomena of the detonation structure. None of this detail 
is apparent at this mesh size in any of the contours. No shear layers were 
observable. Metacomp, maker of CFD++ stress that their Euler solver settings 
employ a behind the scenes dissipation parameter built in the code. We could 
not verify whether this had any effect. 



Figure 36. Temperature contour plots of detonation front using a 100 pm and 

50 pm structured mesh 
50 














Inspecting H 2 mass fraction contours also show scarce signs of unburned 
pockets of hydrogen right behind the front. Based on the literature, this has been 
observed experimentally and has been confirmed computationally. Lee explains 
that local fluctuations of the main shock front during the low velocity phase of the 
cycle can increase induction times resulting in later heat release after the front. 
Many irregular fuels can have significant levels of unburnt reactant. 

C. CHEMICAL KINETICS 

Initial runs used the Kurkov and Burrows chemistry set for hydrogen as an 
initial investigation to see whether a chemistry set developed for supersonic low 
pressure combustion could be utilized in a detonation environment. However, we 
soon discovered this set was not the most ideal set, as diffractions set forth into 
unconfined areas would not completely fail with our smallest structured mesh. 
Granted it is also suspected that this may be a product of an over reactive 
mixture at the current grid resolution. Further investigation was deemed required 
at this point. 

The Hanson and Peterson set was the last iteration, this reduced set was 
exclusively created for high pressure detonation for ram accelerator type 
combustion events. Unfortunately, results were not as promising as hoped, since 
an average was taken for 4 reaction equations, and since they were non- 
Arrhenius and could not be readily inputted into CFD++. These simulations also 
were unsuccessful in producing a failure in the diffraction front. 

Detonation velocity between the various chemistry sets was compared 
and determined for a number of runs based on a finite difference method by 
measuring the front change in position between two time steps, see Table 4 
some the typical values calculated. Values were on the order of 5%-10% 
percent higher than the CEA calculated values. This was expected as the 
nitrogen chemistry was omitted for all Kurkov & Burrows chemistry sets. Our 
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modified Kurkov set included a one step NO equation which was hoped to slow 
down the chemical kinetics. Velocities were slower but the code still was not 
able to resolve a failing diffraction. 


Table 4. Calculated detonation velocity with CEA calculated values 


Run # 

Mesh 

4) 

Reduced 

Numerical Scheme & other 

Avg 


Size 

(mm) 


Chemistry 

Set 

settings 

Detonation 

Velocity 

(m/s) 

6- 

27aug 

0.05 

1 

Kurkov 

Burrows 

Implicit, CFL 1, LHS dissipation 

2010 

7-30 

0.05 

1 

Kurkov 

Implicit, CFL 1, LHS 

2000 

Aug 



Burrows 

dissipation, 1C 1.5Pcj, viscosity 
on walls 


11- 

0.05 

1 

Modified 

Implicit, CFL 1, LHS 

1910 

lOsep 



Kurkov 

Burrows 

dissipation, 1C 2Pcj, added NO 
chemistry 


CEA 


1 


NASA Chemical Equilibrium Vq 

1964 

CEA 


0.8 


NASA Chemical Equilibrium Vq 

1816 


D. DETONATION DIFFRACTION CONTOUR PLOTS 

All simulations exhibited critical diffraction characteristics even though all 
conditions were sub critical. All detonations undergo some partial failure which 
starts almost immediately at the open step but inevitably a re-ignition event 
occurred at some point, reenergizing the entire detonation. Pressure, 
Temperature and OH mass fraction contour plots at one mid time step attest to 
the some of the true complexity of the flow system. Near the main axis of the 
channel the flame front never really appeared to decouple from the precursor 
shock. However, the flame front at the corner completely decouples. Front 
speeds were observed as low as O.SVq near the corner to as high as 1.2Vcj and 
do not ever seem to slow down at the core from the diffraction event. Note the 
vortex at the as seen by the by the streamlines due to rapid lateral expansion. 
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Figure 37. From left to right H 2 , OH mass fractions, pressure (in Pa) 

temperature (in K) 

The following two figures chronicle the detailed re-ignition of the front for 
the standard case with time sequenced pressure and temperature contours. As 
the wave works its way past the corner, a rarefaction wave originates almost 
immediately from the open corner begins to works its way down at an angle of 30 
degrees and works to begin a decoupling of reaction front to precursor shock. By 
frame 6, the rarefaction wave was fully into the reaction zone and sweeps across 
the front starting to flatten it and adding curvature. 

All simulations exhibited a decoupling shock and reaction front from the 
diffraction front followed by a reinitiation event as seen in the above Figures. It 
should be noted that even though reinitiation occurs, the shock never completely 
decoupled. We really only observed localized explosions resulting from 
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rarefaction and transverse wave collisions. This event then triggers a reaction 
across the unreacted zone or shell. Recoupling of the shock and reaction front 
then spreads over the entire front. Before any localized explosion, unburned 
reactants slightly affected by a precursor shock find themselves in a spherical 
zone between a flame front from behind and a transverse wave above. In this 
reinitiation event, it appears the detonation proceeds transversely across the 
chocked reactants in an azimuthally and polar manner of direction completing the 
reaction and reinitiating the entire front, as seen in both the time sequenced 
contours. Figures 38 and 39. 
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Figure 38. Time sequence pressure contour plot with streamlines (in Pa) 
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Figure 39. Time sequenced temperature contour plots (in K) 
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VI. EXPERIMENTAL RESULTS 


A. TESTING OVERVIEW 

Table 5 depicts an overview of all experiments conducted with the split 
flow multi-stepped configured 1:3.5 aspect ratio channel. The tests comprised of 
detonating hydrogen and air at atmospheric conditions for various diffraction area 
ratios. 


Table 5. Experimental test matrix 


Test# 

Test Condition 

Equivalence 
Ratio range 

Diagnostic 

1 

No Step 

0.75-1.5 

Shadowgraph, PIV 

2 

1 Step 1:2.6 area ratio 

0.75-1.5 

Shadowgraph 

3 

1 Step 1:2.3 area ratio 

0.75-1.5 

Shadowgraph 

4 

1 Step 1:2.00 height ratio 

0.75-1.5 

Shadowgraph, 

5 

2 Steps 1:2.3 & 1:2.3 area 
ratio 

0.75-1.5 

Shadowgraph, 

6 

2 Steps 1:2 & 1:2 area ratio 

0.75-1.5 

Shadowgraph, natural 

emmision 


B. PRESSURE DATA 

Low speed pressure measurements were primarily used to calculate 
reactant fuel to air ratio for all tests. With pressure data, equivalence ratio was 
calculated based on the choked mass flow Equation 19, with main air flow using 
a 0.235 inch diameter orifice plate while a 0.157 inch orifice was used for fuel 
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By obtaining a ratio of mass flow rates for both air and hydrogen, the equivalence 
ratio could be determined with reasonable errors of about ±5%. Discharge 
coefficient Cd for both orifice plates was found to be 0.61 using the CRC 
handbook for mechanical engineers. These results were also in good agreement 
with measurements derived PIV flow rates. 

High-speed pressure probe data was acquired using Labview and 
detonation speeds were determined by measuring the inter peak distance with 
time taken into account the probe spacing of 0.241 m. Typical plots obtained can 
be seen in Figure 40. Typical wave speeds ranged from 1,800±50 m/s to 
2,100±50 m/s and agreed very well with wave speed calculated from high speed 
camera imagery and CJ wave speeds calculated from CEA. 



Figure 40. High-speed pressure data 
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C. RESULTS OVERVIEW 

Figures 41 and 42 represents a summary of all tests performed. The 
channel performed best at an area ratio of 2 for a wide range of equivalence 
ratios for both single and double step configurations. It was marginal for the ratio 
of 2.3 and did not transmit for any area ratios of 2.6. We present the data in 
terms of area ratios since we are dealing with a channel that is not perfectly flush 
with the optical section, its support members are a half-inch thick all around so 
there is some diffraction occurring in the z-direction. 
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Figure 41. Single Step diffraction results with Channel Area Ratio 

Detonations transmit into the double step area ratio of 2 for three of the six 
equivalence ratios tested and did not transmit using a double 2.3 ratio, seen in 
Figure 42. 
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Double Step Diffraction Results 
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Figure 42. Double step diffraction results with Channel Area Ratio 


Based on shadowgraph imagery, centerline detonation velocity relative to 

CJ velocities was calculated with respect to distance from the first diffraction 

plane for a number of equivalence ratios and presented on 6 of the following 

Figures. Figures 43, 44 and 45 depict wave velocities just after the first 

diffraction plane at three equivalence ratios ranging from lean to rich with our one 

step configuration for varying height ratios. The area ratio of 2 performed well 

under all conditions as observed to reinitiate well before the second diffraction 

plane (located 3.75 x/h main channel heights away) coming out of the first 

diffraction plane as wave speed were close to CJ velocity. After the second plane 

the wave diffracted into an open space and maintains good strength for another 2 

height distances. Area ratios of 2.3 did not fare as well, as the reinitiation seemed 

to occur at a later distance than 3.75 x/h. Looking at the 0=1 case there was a 

reinitiation that happened, but at that point the wave did not have enough time to 

completely form and gain enough strength to survive once it crossed the second 
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diffraction plane. It seemed that a good step height to length ratio was 2 times 
x/h, so long as the wave survives the initial expansion event. For our area ratio 
case of 2.3, the second step height is 2 inches and so our transition length 
should have been at least 4 x/h, allowing enough distance for a detonation to 
properly reform. Our 2.6 case did not survive the first diffraction as the 
expansion waves were given enough time and distance to penetrate the entire 
front enough precluding any recovery, in all reinitiation cases. 
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Figure 43. Detonation velocity vs. run distance for different area ratios 
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Figure 44. Detonation velocity vs. run distance for different area ratios 
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Figure 45. Detonation velocity vs. run distance for different area ratios 
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Figures 46, 47 and 48 represent the results for both double step setups of 
Area ratio of 2 and 2.3. The double area ratio of 2 successfully transitioned both 
steps with velocities slightly higher in the 3.75-6 x/h zone then the previous case. 
The addition of another wall to created another shock to high impedance material 
interaction creating a Mach stem reinitiation zone due to the shock reflection. It 
should be noted that the range of equivalence ratios for survival was reduced. 
The double area ratio of 2.3 was not as successful. Longer step lengths to 
accommodate this ratio it may have resulted in successful reignition. 
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Figure 46. Detonation velocity vs. run distance for different area ratios 
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Figure 47. Detonation velocity vs. run distance for different area ratios 
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Figure 48. Detonation velocity vs. run distance for different area ratios 
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D. DIFFRACTION PREDICTION MODEL 

A simple diffraction failure model was tested and validated with these 
experimental results. Based on previous results both detonation cell size for 
hydrogen and air with respect to equivalence ratio and detonation velocities 
could be predicted for the various reactants as seen in Table 6. Knowing these 
parameters, our detonation channel's geometric dimensions and how critical our 
channel was for an aspect ratio of 3.5 we could predict the velocity deficit under 
these conditions [24]. We could also calculate either rarefaction angle from 
disturbance to undisturbed detonation a or a reaction zone velocity Vrz to the first 
step from Skews method [25], depending on what data was available. In our 
case we calculated Vrz using Equation 20 and measured a based on the Shultz's 
high speed shadowgraph images of Hydrogen and air and our own shadowgraph 
images which was found to be 15 degrees [8]. Knowing our step distance and 
the distance required for total failure of the front, we could predict the health of 
the front by the time the wave hits the step wall. Degree of criticality in Table 6 
refers to how many detonation cell widths we were across the cross section in a 
particular case with respect to the critical height for a 3.5 aspect ratio channel 
which is 7A. For most cases, we were about 5 times smaller than this so called 
critical length. With this information, we could apply a velocity deficit linked to 
criticality [7]. 
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Figure 49. Schematic of a diffracted shock (From [8]) 


Based on our experiments, it appeared that once the entire front had 
curvature the diffracted chock had no chance of survival. However, if some 
portion of the reaction zone was still coupled to the pre-cursor shock there was a 
chance of survival if a Mach stem developed quickly enough from the wall area 
and had time to progress into the rest of the front. 

Our case with an area ratio of 2.6 failed every time experimentally and 
was predicted in our model. We also find experimentally that the front has 
already been completely disturbed, seen in Table 7. Whereas area ratio of 2 
detonations survived every single step expansion experimentally and we predict 
that this wave was about 80% failed. 
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Table 6. 


Predicted numerical results with a=15°, Acr=7 


<p 

A (in 
cm) 

A's across 

area 

Acr/A's 

across 

(criticality) 

area 

ratio 

(AR) 

criticality 

/AR 

0.8 

2 

1.25 

5.60 

2.00 

2.80 

0.8 

2 

1.25 

5.60 

2.30 

2.43 

0.8 

2 

1.25 

5.60 

2.60 

2.15 

1 

1.5 

1.67 

4.20 

2.00 

2.10 

1 

1.5 

1.67 

4.20 

2.30 

1.83 

1 

1.5 

1.67 

4.20 

2.60 

1.62 

1.15 

1.7 

1.47 

4.76 

2.00 

2.38 

1.15 

1.7 

1.47 

4.76 

2.30 

2.07 

1.15 

1.7 

1.47 

4.76 

2.60 

1.83 


Table 7. Predicted numerical results cent. 


<P 

area 

Vcj 

V deficit 

V 

Time t to 

% 

Exp 


ratio 

(m/s) 

(m/s) 

(m/s) 

wall (s) 

failed 

results 

0.8 

2.00 

1866 

1810.02 

905.01 

4.88E-05 

84.19 

passed 

0.8 

2.30 

1866 

1810.02 

905.01 

5.61 E-05 

96.76 

passed 

0.8 

2.60 

1866 

1810.02 

905.01 

6.35E-05 

109.52 

failed 

1 

2.00 

1966 

1946.34 

973.17 

4.54E-05 

84.19 

passed 

1 

2.30 

1966 

1946.34 

973.17 

5.22E-05 

96.76 

passed 

1 

2.60 

1966 

1946.34 

973.17 

5.91 E-05 

109.52 

failed 

1.15 

2.00 

2015 

1974.7 

987.35 

4.48E-05 

84.19 

passed 

1.15 

2.30 

2015 

1974.7 

987.35 

5.15E-05 

96.76 

passed 

1.15 

2.60 

2015 

1974.7 

987.35 

5.82E-05 

109.52 

failed 


E. TEST CONDITION 1: COLD FLOW AND WIDE OPEN DIFFRACTION 

Test condition 1 included all experiments that would help understand all 
the conditions the detonation front would be exposed to in the simplest 
configuration, i.e., diffraction into a wide opening without the stepped sections in 
the optical section. Tests included basic PIV tests of the cold flow and 
diffractions into an open channel as well shadowgraph shots. 
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Cold flow tests were conducted to better understand the pre detonation 
flow conditions and seeding quality of the flow. Figure 49 was a typical test run 
with main air flow mass flow rate of 0.017 kg/s depicting camera bit count 
(intensity level) of the seed. The laser sheet illuminates an 80 mm by 80 mm 
zone in the flow immediately following the diffraction plane, flow moves from left 
to right. Note that these counts were slightly high as the bulk of the seed was 
biased towards lower main detonation channel. Although seed density was 
higher in this region than required at these low speed fills, it was found that such 
levels were required for sufficient seed density required for good vector 
correlation for wave speeds in the order of 2,000 m/s as such waves would tend 
to disperse the seed after the front, as particles can be accelerated to as much 
as 1,000 m/s just after the front applying significant seed density variations. 
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Figure 50. Seed bit counts with low speed cold flow run 


The postprocessed vector plot for Figure 50 can be seen at Figure 51. 
Flow on average was very turbulent and fluctuated significantly during the filling 
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process. On average, velocities in the top channel were higher by a factor of 2 
compared to flow in the lower channel. A shear layer of lower velocity can be 
found in between both channels. 
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Figure 51. Cold Flow PIV velocity vector plot 


The following Figures depict PIV vector plots of one snapshot of the 
detonation front without any steps with a 0=1. The first Figure is frame two of 
the two frame camera showing the good seed density in and around the 
diffracting detonation front. A pressure of 5 psig was applied to the seeding 
system for reasonable seeding density. Inter laser time was calculated in Davis, 
Lavision's in house PIV code based on estimated 800 m/s particle velocity behind 
the front, resulting in a 710 ns time deference between frames dt. The trigger 
delay was set to 230 ps and triggering was off of the first Kistler pressure probe 
with a laser power of 50% of max. 
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velocity m/s 


















Figure 52 was the post processed vector plot of the flow field. We 
observed particle speeds in the 800 m/s range near the core dropping to as low 
as 400 m/s if we follow around the perimeter. Unfortunately, with such a small dt 
we missed capturing lower end velocities and similar test runs would be needed 
with various dt's to get the whole event. Additionally, because we used a low 
frame rate PIV speed system to capture detonations at various diffraction times 
many runs would have been required. 
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Figure 52. Seed bit counts with high speed settings 


The postprocessed plot was quite noisy around the front. It is suspected 
that fine tuning camera focus and seed density would have improved the 
proportion of noise and we would have obtained better correlation with calculated 
vectors. 
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Figure 53. Post processed vector PIV vector plot of a diffracting detonation 

front 


Shadowgraph sequences were also taken without any step sections this 
can be seen in Figure 54. All images were shot at 150,000 frames per second 
with a total of 17 frames with the Ultra 17 high-speed camera system. 
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F. 


TEST CONDITION 2: SINGLE STEP AREA RATIO 2 WITH 
SHADOWGRAPH 



Figure 55. One step, 0=1,2.00 A 2 /A 1 



Figure 56. Frame 4 From above, x/h 1-3.5 
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G. 


TEST CONDITION 3: SINGLE STEP AREA RATIO 2.3 WITH 
SHADOWGRAPH 



Figure 57. One step, 0=1,2.30 A 2 /A 1 



Figure 58. Frame 8 from above, x/h 2-4.5 
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H. 


TEST CONDITION 4: SINGLE STEP AREA RATIO 2.6 WITH 
SHADOWGRAPH 



Figure 59. One step, 0=1,2.6 A 2 /A 1 



Figure 60. Frame 9 from above, x/h 2 - 4.5 
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I. 


TEST CONDITION 5: DOUBLE STEP AREA RATIO 2.00 WITH 
NATURAL EMISSION 


The following two images were taken with high-speed camera looking 
directly at the optical section without laser illumination. A reinitiation event can 
be observed moving azimuthally across what was once a decoupled reaction 
zone. In Figure 62, the front is still planar at a distance of 5.5 x/h. 



Figure 61. Natural emission, 0=1,2.00 A 2 /A 1 , x/h 1 - 6 



Figure 62. Natural emission, 0=1,2.00 A 2 /A 1 , x/h 1 - 6 
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J. 


TEST CONDITION 6: DOUBLE STEP AREA RATIO 2.3 WITH 
SHADOWGRAPH 



Figure 63. Two step, 0=1,2 A 2 /A 1 



Figure 64. Frame 12 from above 
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VII. FUTURE WORK 


Future work could encompass both computational and experimental 
aspects. Computationally, the next step would be to rerun these results with a 
more refined structured mesh and to develop a better chemistry set for 
detonations. 

Experimentally, longer step sections could be manufactured to ensure 
proper full transition during diffractions events. Also, different fuels could be 
used to test the predictions model as well as a performance comparison to 
hydrogen. Additional PIV work is required once Lavision rebuilds the lab's 
existing camera. With the new optical section more diffraction experiments could 
be conducted using CH* and PLIF. 
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VIII. CONCLUSIONS 


A new optically-accessible test section was designed and built to image 
the detonation diffraction of hydrogen and air mixtures across multi-step 
configurations. Shadowgraph and PIV imaging techniques were used to 
determine shockwave profiles, velocities, and re-ignition mechanisms throughout 
the diffraction process. Computational models were investigated but failed to 
produce any useful information due to the artificial high reactivity of the mixtures 
for the various reduced chemistry kinetics and also inadequate mesh resolution 
utilized. 

Experiments revealed that hydrogen/air detonations could survive single 
and double diffraction across an area ratio of 2.0 for an equivalence ratio range 
of 0.8 to 1.14. Successful detonation diffraction was observed for area ratios of 
2.3 but only for the primary diffraction event and only for equivalence ratios near 
stoichiometric for the given step length. Increasing the step length would 
improve the likelihood of subsequent diffraction events by allowing the detonation 
to reform before encountering another diffraction condition. No successful 
detonation diffraction reinitiations were observed for area ratios of 2.6. 
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APPENDIX A. NASA CHEMICAL EQUILIBRIUM OUTPUT 


NASA-GLENN CHEMICAL EQUILIBRIUM PROGRAM CEA2, FEBRUARY 5, 2004 

BY BONNIE MCBRIDE AND SANFORD GORDON 
REFS: NASA RP-1311, PART I, 1594 AND NASA RP-1311, PART II, 1996 


-dr lA- rJl- -dr -A; it :A- -dr ^ ^ -A: ^ :A- it ^ ^ -A: it ^ Jr Jc Jc Jr Jr Ji; Jc ■i' -Jr Jr Jr Jc Jc -A -Js Jr Ji; Jc, -Jr -Jr Jr Jc Jc -Jr Jr 

Jc ik ik Jr Jr Je Je ■i' 


output massf 

prob ca3e=000u2 142 det 

phi=l 

p [atm) =1 
output short 
output trace=le-iO 


reac 

oxid A.ir t, k:= 298 wt%= 100. 
fuel H2 t, k:= 2 98 wt%= 100. 
end 

DETONA.TION PROPERTIES OF AN IDEAL REACTING GAS 

CA.SE = 00002142 


TEMP 


REACTANT 


K 

OXIDANT Air 

298.000 

FUEL H2 

298.000 


WT FRACTION 
(SEE NOTE) 

1.0000000 
1 .0000000 


ENERGY 
KJ/KG-MOL 
-129.855 
-4 . 325 


0/F= 34.29623 %FUEL= 2.833164 R,EQ.RATIO= 1.000000 PHI, EQ . RA_TIO = 

1 .000000 


UNEURNED GAS 


PI, BAR 1.0132 
Tl, K 290.00 
HI, K.J/KG -4.36 
Ml, (1/n) 21.008 
GAMMAl 1.4015 
SON VELl,M/SEC 406.6 


BURNED GAS 
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T, K 

2943.63 

RHO, KG/CU M 

1.5500 0 

H, KJ/KG 

1333.33 

U, KJ/KG 

313.90 

G, KJ/KG 

-29833.7 

S, KJ/[KG){K} 

10 . 5080 

M, (i/n) 

24.008 

[dLV/dLP)t 

-1.00954 

CdLV/dLT)p 

1 . 2062 

Cp, KJ/[KG)(K) 

3.3542 

GAMKAs 

1 . 1637 

SON VEL,M/3EG 

1009.2 

DETONATION PARAMETERS 

P/Pl 

15.595 

T/Tl 

9.078 

M/Ml 

1 . 1428 

RHO/RHOl 

1 . 0042 

DET MACH NUMBER 

4.0335 

DET VEL,M/SEC 

1965.1 


MASS FRACTIONS 


*ALr 

1.2550-2 

*CO 

1.3093-4 

*C02 

2.6525-4 

COOH 

2 . 6149-9 


2 . 4760-4 

HCO 

1.0225-9 

HNCO 

2.395-10 

HNO 

3 . 1235-6 

HN02 

8 . 9839-7 

HN03 

1.020-10 

HO 2 

1 . 5250-5 

*H2 

2.6254-3 

HCOOH 

3 . 735-10 

H20 

2 . 2035-1 

H202 

2 . 3000-6 

*N 

1 . 1146-6 

^NH 

4.5040-7 

NH2 

2 . 7000-7 

NH3 

3 . 4969-7 

NH20H 

4.217-10 

*NO 

9.4235-3 

NO 2 

6 . 6339-6 

*N2 

7 . 2930-1 

N20 

2 . 7711-6 

*o 

1 . 3571-3 

*OH 

1.3522-2 

*02 

1 . 0117-2 

03 

6.2070-9 
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APPENDIX B. CAD DRAWINGS 


The following set of Figures were all generated with the Solidworks 2010 
mechanical modeling software package. Optical section and the new head end 
inlet flange were made from Type 303 Stainless Steel in order to meet the rigidity 
and strength requirements for the engine test section while also being the most 
readily mechineable? for an austenitic stainless steel. All internal components 
were made from 6061 Al. All units are in inches and degrees. Tolerances for 
linear dimensions are ±0.01" and radial dimensions are ±0.05°. 
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Inlet Manifold 



Optical Assembly 


Figure 65. Overall detonation channel assembly 
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BK7 glass 


Back Plate 


Glass Holder 


Top plate 
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Figure 66. Exploded view of Optical assembly 
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Figure 67. Back plate (1 of 2) 


88 






























































8 .^ 




Figure 68. Back plate (2 of 2) 
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Figure 69. Outer wall optical left and Right (1 of 4) 
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DETAIL D 
SCALE 1 : 1 
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Figure 71. Outer Wall optical Left and Right (3 of 4) 

92 














































































































































































































Figure 72. Outer Wall optical Left and Right (4 of 4) 
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Figure 76. Glass Holder optical (2 of 2) 
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Note 1 




Note 1: optical polished surface or 
both front and back 

Note 2: all sharp edges and corner: 
removed 


Figure 78. BK7 glass for main optical section 











Figure 79. Bottom plate assembly for PIV laser sheet 
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Figure 80. Bottom plate with BK7 glass slot PIV laser sheet (1 of 2) 
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SECTION D-D 
SCALE 1 :2 
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Figure 84. Aluminum Insert 
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Figure 88. Insert top (1 of 2) 
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DETAIL C 
SCALE 1 :4 


Figure 89. Insert top (2 of 2) 










Figure 90. Insert bottom (1 of 2) 
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DETAIL 6 

SCALE 1 :2 


Figure 91. Insert bottom (2 of 2) 
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Figure 92. Insert left (1 of 2) (mirrored for right side) 
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Figure 93. Insert left (2 of 2) (mirrored for right side) 




















Figure 94. Insert front 
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Figure 95. 


Insert front adapter 
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SECTION D-D 
SCALE 1 : 2 


Figure 96. Insert ramp holder (1 of 2) 
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Figure 97. Insert ramp holder (1 of 2) 


118 
















00.11 ^ 0.10 




Figure 98. Small ramps 
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Figure 99. Step insert sample (steps made for 1:1.25, 1:1.5 and 1:1.75) 
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Figure 102. Inlet manifold (1 of 2) 
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APPENDIX C. STANDARD OPERATING PROCEDURE 


Test Cell #1 Standard Operating Procedures (SOP) 

Combustor Start Up 
(Modification Date 4 Nov 10) 

INITIAL SAFETY CHECKS 


1. Notify all lab personnel of intention to make test cell 1 live. 

2. Turn ON control console 

3. Turn ON warning lights 

4. Cell #1 EMERGENCY SHUTDOWN BUTTON (Control Room)- VERIFY 
IN 

INITIAL PREPARATIONS 

1. Notify the Golf Course (x2167) (Only required if Hot Fire Test is 
conducted) 

2. Open Test Cell Door 

3. Igniter Control (Test Cell)-VERIFY OFF (Red Button Out) 

4. PXI-1OOOB Rack (Test Cell)-VERIFY ON 

5. Shop Air-VERIFY > 100 PSI 

6. Shop Air Valve (Test Cell)-VERIFY OPEN 

7. 115 VAC Control/Cell #1 Switch (Control Room)-ON 

8. 28VDC Power Supply/Cell #1 Switch (Control Room)-ON 

9. Open LABVIEW file "TCI ControllerMike" 

TESTING SET-UP 

1. Inspect optical section and mirrors for cleanliness 

2. Commence Argon-Ion laser start-up procedure (see Arg-lon SOP) 

3. Commence PIV NG-Yag double pulse laser start up (PIV SOP) 

4. PIV software controller - ON 

5. Kistler Amplifier Power-ENSURE OFF 

6. Kistler Leads-CONNECT 

7. Exhaust Tube-VERIFY PROPER POSITION 

8. Notify all personnel that gasses and TESCOM will be enabled. 

9. Test Cell #2 and #3 Node 1 Air Isolation Valve (Test Cell #2)-VERIFY 
CLOSED 

**NOTE: This valve maybe left open only if Test Ceil #2 or #3 is configured 

for active testing** 
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10. TRANSDUCER and TESCOM Power Switch (Test Cell #2)-ON 

11. Set 0 (Zero) pressures on ER3000 (Control Room) for the following: 

a. Node 3 (Hydrogen) 

b. Node 1 (Main Air) 

c. Node4 (Seeder Air) 

12. Main HP Air Jamesbury Valve (Outside Test Cell)-OPEN SLOWLY 

13. Power Strip (above PXI-1000B)-VERIFY ON 

14. Igniter Control Light 5V logic (Red LED CRYDOM 2nd Column, 2nd down 
in PX-1000B Rack)-VERIFY OUT 

**DANGER: IF RED LIGHT IS ENERGIZED, MUST CLEAR USING LABVIEW 
BEFORE CONTINUING TO PREVENT PREMATURE IGNITION** 

15. Igniter Control (Test Cell)-PUSH RED BUTTON IN 

16. Igniter Control Startup Diagnostics-OBSERVE COMPLETION OF 
DIAGNOSTICS (Verify energy level setting reads 2.03 J) 

17. Main HP Air Isolation Valves (2) (Located in Test Cell)-OPEN SLOWLY 

**DANGER: OPEN VALVES SLOWLY TO PREVENT RAPID 
PRESSURIZATION OF DOWNSTREAM LINES** 

18. Node 3 (Hydrogen) Shop Air Valve (Above Bottle in Bottle Room)-OPEN 

19. Hydrogen Bottle Isolation Valve (On Bottle)-OPEN SLOWLY 

** VERIFY ADEQUATE ETHYLENE PRESSURE FOR TESTING ON 
DOWNSTREAM GAGE LOCATED IN BOTTLE ROOM** 

PRE-TESTING 

1. SA-5 Camera-TURN ON AND REMOVE COVER 

2. Start SA-5 Software 

3. Verify Image on SA-5 Software 

4. PIV Camera - TURN ON AND REMOVE COVER 

5. Start PIV software 

6. Verify Desired Trigger Types and Valve Durations 

7. Determine Desired Fuel and Air Pressures (Mass Flow Choke 

Calibrations.xls) 

8. Set Required Pressures on ER3000 

a. Node 3 (Hydrogen) 

b. Node 1 (Main Air) 

c. Node 4 (Seeder Air) 

9. Ready Digital Camera 

10. Switch Monitor to Test Cell #1 
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TESTING 


1. Clear All Test Cells and Verify with Head Count 

2. Flashing Red Light and Siren-ENERGIZE 

3. Verify Golf Course is CLEAR 

4. Digital Camera RECORD 

5. Camera-ARM 

6. In LABVIEW Enable Facility Button-ON 

7. Test Cell #1 Emergency Shutdown Button-TURN CLOCKWISE 

8. In LABVIEW Start Button-CLICK TO START 

WHEN TESTING COMPLETE 

9. Set Node 1 Pressure to 0 (Zero) 

10. In LABVIEW Turn Off Button-CLICK TO SECURE 

11. In LABVIEW Enable Facility Button-VERIFY OFF 

12. Test Cell #1 Emergency Shutdown Button-PUSH IN 
13.Siren-OFF 

14. Digital Camera-Stop/Pause 

15. Save SA-5 Image 

16. Save PIV data 

**NOTE: If Further Testing is Required, re-perform steps 1-15 of the 

Testing Section** 

POST-TESTING 

1. Close fuel isolation valves, rerun labview code to cycle fuel valve 

2. Igniter Control (Test Cell)-PUSH RED BUTTON OUT 

3. Flow additional Main Air to purge any remaining fuel 

4. Set ER300 Node 3 (Ethylene) to 0 (Zero) 

5. Verify pressures on ER3000 are set to 0 (Zero) on the following: 

a. Node 3 (Hydrogen) 

b. Node 1 (Main Air) 

c. Node 4 (Seeder Air) 

6. Main HP Air Isolation Valves (2) (Located in Test Cell)-CLOSE 

7. Main HP Air Jamesbury Valve (Outside Test Cell)-CLOSE 

**NOTE: If Further Testing will be accomplished with a different Ramp 
obstacie configuration, return to the TESTiNG SET-UP Section, if not 

continue to step 7.** 

8. SA-5 Camera-TURN OFF AND INSTALL COVER 

9. PIV Camera - TURN OFF AND INSTALL COVER 
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10. TRANSDUCER and TESCOM Power Switch (Test Cell #2)-OFF 

11. Close Test Cell Door 

12. Node 3 (Hydrogen) Shop Air Valve (Above Bottle in Bottle Room)-CLOSE 

13. Secure Bottle Room 

14. Exit out of SA-5 and PIV software 

15. EXIT out of LABVIEW 

10.28VDC Power Supply/Cell #1 Switch (Control Room)-OFF 
11.115 VAC Control/Cell #1 Switch (Control Room)-OFF 
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